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The propagation of Langmuir waves in plasmas is known to be sensitive to density fluctuations. Such 
fluctuations may lead to the coexistence of wave pairs that have almost opposite wave-numbers in the 

O i 1 vicinity of their reflection points. Using high frequency electric field measurements from the WIND 

satellite, we determine for the first time the wavelength of intense Langmuir wave packets that are gener- 
ated upstream of the Earth's electron foreshock by energetic electron beams. Surprisingly, the wavelength 
is found to be 2 to 3 times larger than the value expected from standard theory. These values are consistent 

<**> ^ with the presence of strong inhomogeneities in the solar wind plasma rather than with the effect of weak 

&) ■ beam instabilities. 
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^y 1 Introduction 

The theory of beam-plasma interactions in homogeneous plasmas is one of the pillars of plasma physics, 

and yet it fails to properly describe some important physical phenomena observed, such as the generation 

of type III solar radio bursts or waves registered in the electron foreshock region, upstream of the Earth's 

qv ■ bow shock (Ginzburg, 1970; Robinson, 1997). The major reason for this discrepancy between observations 

£f) | and the theory of quasilinear, weak or strong turbulence is supposed to be related to the existence of strong 

density fluctuations that are not associated with the beams. When these fluctuations are large enough, 

Langmuir waves can be trapped or reflected and in both cases the wave-vector k\ of the reflected wave 

f*»» ■ becomes approximately equal and opposite to that of the incident wave fco — —k\. Such strong density 

fluctuations are not resolved by existing particle detectors. However, they provide a useful means for 

inferring the wavelength, which is difficult to measure directly in space plasmas. 

The WAVES instrument (Bougeret et al., 1995) onboard the WIND satellite routinely measures Langmuir 
waves at high time resolution (up to 120000 samples per second) in the solar wind. When two large 
amplitude and counter-propagating waves coexist, a monochromatic wave packet with a slowly varying 
envelope is produced. Some typical examples are shown in Fig. 1. Here, for the first time, we show how 
the wavelength of these waves can be directly estimated from a polarization analysis of two components of 
the electric field. Interestingly, the measured value disagrees with estimates derived from standard beam- 
plasma theory in homogeneous plasmas but is consistent with estimates derived for strongly inhomogeneous 
plasmas. 



2 The effect of density fluctuations 

Density fluctuations in the solar wind plasma have been studied for many years (Neugebauer, 1975; 
Celnikier et al., 1987) and exhibit a power law over many decades. The smallest scales are difficult to mea- 
sure by lack of time resolution and counting statistics of standard particle instruments. Recently, Ergun et al. 
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Figure 1: Five typical coherent Langmuir wave packets observed by WIND. Only the E x electric field 
component is shown. The bottom plot shows a marginal case in which the wave-packet is barely coherent. 

(2008) observed the interaction of density fluctuations with Langmuir waves, and described them as eigen- 
modes trapped inside density depletions. In such strongly inhomogeneous plasmas, wave propagation is 
strongly affected by density fluctuations. As a consequence, the dynamics of the instability departs quite 
substantially from that observed in homogeneous plasmas (Nishikawa and Riutov, 1976; Muschietti et al., 
1985; Krasnoselskikh et al., 2007). Such density inhomogeneities can cause high frequency electrostatic 
waves to be converted into electromagnetic waves (Bale et al., 1998; Kellogg et al., 1999). They also af- 
fect the statistics of wave amplitudes (Cairns et al., 2000; Krasnoselskikh et al., 2007) and the presence of 
sufficiently large inhomogeneities eventually leads to wave reflection (Ginzburg, 1970; Budden, 1988). 

Here, we focus on the dynamics of Langmuir waves that are generated near the electron foreshock edge 
region by electron beams with a typical energy of 1-3 keV (Filbert and Kellogg, 1979). These waves have 
a narrowband spectrum and their characteristic frequency is close to the electron plasma frequency. The 
presence of relatively large but slowly varying (with respect to the wavelength) density fluctuations, with 
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5n/n > -k X 



(1) 



can substantially alter the beam-plasma instability by changing the wave-vector k along the trajectory of 
the beam. Here, k is the wave-vector of the plasma wave, tlq the background plasma density and Ad is the 
Debye length. The wave trajectory in the inhomogeneous plasma can be described by the equations 
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that results in a fc-vector change with the density that is described by 
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The observations we present hereafter clearly show the simultaneous presence of two oppositely propa- 
gating waves with fc-vectors close to each other, so that fci ra — fc . Two scenarios can lead to such a 
coexistence of primary and oppositely directed secondary waves. The first one is the scattering off of 
density fluctuations, and the second one is three-wave decay. 

Let us consider the first scenario. Without loss of generality, we assume here that the density increases 
along the x direction. In this case, the k x component of the wave-vector and the wave group velocity both 



decrease along x. The wave amplitude then increases even in the absence of growth or damping. This effect 
can be illustrated in the WKB approximation, where the solution of the wave equation can be written in the 
form 



E ~ EqJ cxp I -jut + j fcx(x) dx I , (5) 

and the amplitude increases when k x (x) decreases. Without damping or growth, the component of the 
energy flux F along the density gradient can be supposed to be constant. We then have F ~ V gi E 2 ~ 
k^E 2 = const and so the wave amplitude E scales like k x ' . In the WKB approximation, reflection 
processes are supposed to be similar to mirroring, with part of the energy being reflected and part damped. 
Two large amplitude waves, incident and backscattered, can thus be observed simultaneously in the vicinity 
of the reflection point (Bale et al., 2000). It is worth mentioning that the WKB approximation is valid close 
to, but not at, the turning point where k x goes to zero. Immediately near this point the wave amplitude 
resembles solutions of the Schrodinger equation with a linearly varying potential expressed in terms of 
Airy functions and it remains finite at the reflection point (see Sect. 17 in Ginzburg, 1970). 

The second scenario is a change in the characteristics of the wave-wave interaction by a variation of the 
wave-vector. The growth of primary waves generated by beam-plasma interactions can be saturated by the 
decay process Iq — > l\ + s, in which a primary Langmuir wave ^o decays into a secondary Langmuir wave 
li and an ion-sound wave s. The threshold condition for the decay instability is 

£ ° g ° > i*i*. (6) 

2n Q k B T w s w L 

Here, E$ is the amplitude of the primary wave, T is the electron temperature, i/g and is^ are respectively the 
damping rates of the secondary ion sound and Langmuir waves, and u>s and wl are their pulsations. When 
the primary wave propagates into a denser plasma, the decrease in k x leads to a decrease of the wave-vector 
of the secondary waves. The phase velocity of the secondary Langmuir wave then grows and the damping 
rate drops, finally initiating a decay instability. It is well known from linear theory that the maximum 
increment of decay instability then corresponds to almost exact back scatter with a small correction to 
comply with momentum and energy conservation laws. As a consequence, in this case the wave-vector 
k\ of the secondary Langmuir wave becomes also approximately equal and opposite to that of the primary 
wave, i.e. k\ ~ — ko. 

These oppositely-directed waves provide a unique opportunity for measuring their wavelength without re- 
quiring multipoint measurements. 



3 Expressions for evaluating the wavelength 

We assume that the frequency difference 5f = /o — /i (where /o and f\ are respectively the frequencies 
of the incident and reflected waves) is mainly caused by the Doppler shift rather than by the frequency 
difference associated with the differing dispersions of the two waves. Indeed, this frequency difference 
significantly exceeds that of the ion-sound wave that might be generated by the decay process. In the 
spacecraft reference frame, the wave frequencies are 

i 2 \ 2 \ , fa' ^SW 



fi = f P II + 2 fc i A d ) + -^f 1 - l = ^ C7) 

where V$yj is the velocity of the plasma with respect to satellite. We assume that 

^-^ » UfM, (8) 

uj p 2 

and will check this a posteriori. Theoretical predictions and observations both lend support to the alignment 
of the wave-vector of the primary (and consequently secondary) waves along the magnetic field. Since the 
two wave-vectors have almost equal magnitudes, we find 




where we suppose that both wave-vectors are approximately co-aligned with the magnetic field. From this 
expression, we finally obtain the wavelength A = 2tt /ko. Let us now show how these different quantities 
are estimated. 
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Figure 2: Analysis procedure of a typical wave packet: (a) two components of the electric field in GSE 
coordinates, (b) identification of the two peaks in the power spectral density (continuous curves are gaus- 
sian fits), (c) verification that the amplitudes and the phases (not shown) indeed vary slowly in time, 
(d) hodogram to compare the polarisation with the orientation of the magnetic field. The direction of 
the solar wind velocity varies and is not necessarily aligned with the magnetic field. 

4 Determining the wavelength in practice 



In practice, to measure the wavelength of the wave packets, we first need to know their direction of propaga- 
tion, which is a prerequisite for estimating the wave-vector k. We first convert the electric field components 
into Geocentric Solar Ecliptic (GSE) coordinates, hence defining components E x and E y , see Fig. 2a. Next, 
a continuous Fourier transform is applied in a small interval around the plasma frequency. Except for some 
distorted wave-packets whose modulation width T is only a few oscillation periods, the power spectral den- 
sity usually exhibits two distinct peaks, see Fig. 2b. The relative width A/// (where A/ is the full width at 
half maximum) of these peaks is always well below 1 %, confirming the stability of the frequency over the 
interval. Each component of the electric field can then be modeled as a sum of two complex exponentials. 
A least-squares regression of the following model is performed 

Ei(t) = A ,i(£)Re[exp(j27r/ £+.7<Kx(*))] 
+ A lti (t) Re [cxp (j2irfit + j> M (t))] 
+ Q(t) + e,-(t) 

where i — x or y and e(t) is the residual error to be minimized. The amplitude A(t) and the phase <j>(t) 
vary slowly with respect to T. Possible offsets are taken care of by a small and slowly varying constant 
C(t). We use a sliding gaussian window of width W to estimate these parameters at different parts of the 
wave-packet. Obviously, we must have l//o <C W <C T. Taking W = 10//o was found to be a good 
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Figure 3: Example of electric field spectrum from the WAVES instrument (top plot) and electron flux from 
the Three-Dimensional Plasma and Energetic Particle Investigation (Lin et al., 1995) on 17 April 1996. 



compromise. We checked that the amplitudes Ao(t) and Ai(t) show no beating like the original wave 
packets do (Fig. 2c) and that the phases <f>o(t) and <fii(t) indeed vary slowly with respect to the wave period. 

Once the amplitudes and the phases are known, we can separate the components of the electric field cor- 
responding to the primary and secondary waves. A plot in polar coordinates (see Fig. 2d) reveals their 
polarisation in the xy -plane. Also shown are the projection on that plane of the magnetic field B and the so- 
lar wind velocity Vsw- For electrostatic waves, V-E = 0^-k//E and so the direction of the wave-vector 
k should be given by the semi-major axis of the polarisation ellipse. 

Figure 2d shows two remarkable results, which are systematically found for all the coherent wave-packets 
we analyzed. First, both waves are almost linearly polarized and their polarisation ellipses have the same 
orientation. We conclude that their wave-vectors are necessarily parallel or antiparallel, as expected. Sec- 
ondly, the angle between the semi-major axis of the polarisation ellipse and the ambient magnetic field is 
always small, rarely exceeding 20dcg. The absence of the third E z component in principle prevents us from 
concluding that k // B. For all the events we analyzed, however, we always found the polarisation ellipse to 
be parallel to the magnetic field, no matter how the latter is oriented. We therefore conclude that both waves 
must necessarily propagate along or very close to the magnetic field direction. Momentum conservation 
finally implies fco ~ — k\. 

For the example shown in Fig. 2, we measure Sf = 330 Hz, f p = 20.815kHz and so from Eq. (9) we 
finally obtain the wavelength A = 2.57km. We check that the waves are indeed weakly dispersive, since 
from Eq. (8) we have 

k • Vsw „ nin _ 3 2 2 



= 0.017 > ^fcV> = 0.0031 



(10) 



Similar values are obtained for the other coherent wave packets shown in Fig. 1, except for the last one, 
for which /o and /i cannot be properly identified. We analyzed 644 wave-packets that way and found the 
wavelength to be on average A = 2.1 km, with a standard deviation of 0.5 km; the typical Debye length in 
the solar wind at 1 AU is 10 m. 



5 Discussion and conclusion 

The coherent wave packets we observe always coincide with an enhancement in the energetic electron flux. 
The energy of these electrons ranges from several hundreds of eV up to several keV. Figure 3 shows the en- 
ergetic electron flux measured by WIND during a typical burst of Langmuir wave activity. Fitzenreiter et al. 
(1996) have shown that these wave packets are associated with bump on tail features in the electron distri- 
bution. Their wavelength then should be A = 2irVb/Ld p , where Vb is the beam velocity. For waves with a 
frequency of 20 kHz that are generated by 1 keV electron beams, we find A = 935 m. This value is about 
2.3 times smaller than the one we measure. The two values cannot be reconciled by assuming that the wave 
activity is triggered by higher energy energy electrons, since unrealistic energies in excess of 5 keV would 
be required. 

The wavelength we observe, however, is compatible with that found for wave propagation in inhomoge- 
neous plasma. As it was pointed out by Kellogg et al. (1999) from the analysis of the density fluctuation 
spectra, variations in the /c-vector, as caused by density fluctuations, can exceed the vector magnitude. 
Therefore, the discrepancy we observe between the predicted and observed wavelengths is most probably 
due to a change of the wave-vector when the waves propagate into a higher density region. Indeed, for the 
solar wind conditions we consider, the variation of the wave dispersion can be as large as 

Suj/uj p = (3/2)k 2 \l ra 0.3%. (11) 

The associated level of density variations is 

Sn/n « 2Suj/uj p « 0.6%. (12) 

Such levels are compatible with estimates made at the electron foreshock (Neugebauer, 1975; 
Celnikier et al., 1987; Kellogg et al., 1999) and are large enough to offer an explanation for the relatively 
large wavelengths we observe. 

Our conclusion is further supported by recent direct observations by Ergun et al. (2008) and 
Krasnoselskikh et al. (2007) of unusually high levels of density fluctuations in the solar wind. The con- 
ventional theory of beam-plasma interaction that has been developed for homogeneous plasmas does not 
apply to such cases. Indeed, our results show in an unambiguous way that the observed wavelengths ex- 
ceed their theoretical prediction. The theory should therefore include the effect of large amplitude random 
density fluctuations from the very beginning. 
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